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Abstract 

The finite temperature phase diagram of a U(l) Higgs- Yukawa model at a finite 
value of the scalar self coupling A is investigated by means of a large- Nf calculation and 
numerical simulations. The phase diagram is similar to the one at zero temperature and 
shows a ferromagnetic, two symmetric and an antiferromagnetic phase. However, the 
phase transition lines are shifted to larger values of the Yukawa coupling demonstrating 
the occurence of the finite temperature symmetry restoration. 



1 Introduction 



The finite temperature electroweak phase transition plays an important role in scenarios of 
the early universe which assume that at a high enough temperature the universe was in a hot 
symmetric state 0. During its expansion the universe cooled down and at a critical tempera- 
ture passed through the electroweak phase transition breaking the symmetry spontaneously. 
That the symmetry should be restored at high enough temperatures was already predicted 
in the work of Kirzhnits and Linde . A subsequent analysis |3|, f|] within the framework of 
large N and perturbative approximations confirmed the symmetry restoration picture. To 
test the validity of these proximations, it would be useful to investigate the phenomenon 
of symmetry restoration also non-perturbatively. Moreover, questions like the order of the 
phase transition or the value of the critical temperature T c can not be answered reliably in 
perturbation theory alone. Nonperturbative lattice simulations may provide more insight 
into the nature of the electroweak phase transition. However, the inclusion of all degrees 
of freedom is too demanding with present day computers and one has to restrict oneself to 
subsystems of the standard model like gauge-Higgs systems or pure scalar theories. In 
particular in the 0(4) symmetric 4 theory the critical temperature can be estimated assum- 
ing that the gauge and fermionic degrees of freedom can be neglected. As the gauge coupling 
is rather weak this seems to be a reasonable approximation and one finds T c « 350GeV [0. 

In this letter I want to include the fermions by studying a U(l) Higgs- Yukawa system in 
the large fermion number (Nf) approximation in combination with lattice simulations. Higgs- 
Yukawa models on the lattice have been the object of numerous analytical and numerical 
studies in the last years 0. They are expected to provide insight into non-perturbative 
properties of the Standard model. They revealed, for example, a surprisingly complex phase 
diagram. For a finite value of the scalar self coupling A four different phases were found at 
zero temperature || |9|. These phases are a ferromagnetic phase (FM), where the symmetry 
is spontaneously broken, two symmetric phases (SYM), one at weak and the other at strong 
Yukawa coupling and an antiferromagnetic (AFM) phase where the staggered magnetization 
is non-zero with two regions, separated by a first order phase transition line. For an infinite 
scalar self coupling A = oo even an additional phase was found, a ferrimagnetic (FI) phase 
where the magnetization as well as the staggered magnetization is non-zero [|I~C|j . 

At small scalar self coupling A the system can be analyzed by analytical means like large- 
Nf, perturbative or mean- field approximations. The non-perturbative simulation results are 
in good agreement with these analytical calculations || . In particular it could be confirmed 
that the phase transition from the ferromagnetic to the symmetric phase at weak Yukawa 
coupling, which is relevant for the standard model, is second order with the critical behaviour 
of the Gaussian fixed point. This means the triviality of the Higgs Yukawa models. Of course, 
the model is plaqued by the appearance of doubler fermions so that one can not expect 
quantitative results from these studies. Nevertheless, Higgs- Yukawa models are expected to 
describe qualitative features correctly like the order of the phase transitions and the question 
of symmetry restoration which I want to study here. 
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2 The model and the techniques 



2.1 Zero temperature 

The lattice action for the U(l) chiral invariant Higgs- Yukawa model at zero temperature is 
defined as 

S = S f + S H . (1) 
The fermion part of the action Sf is given by 

Sf = Y,$i(xW(x,z)i/; i (z) , i = l,2,...,N f /2 . (2) 

In eqn. (0) the fermion matrix may be written as 

M(x, z) = 7m [$x+h,z - ^x-n,z] + V [0i {x) + ij5<t>2 (x)] S XiZ , (3) 

where 7 M , 75 are the Hermitian Dirac matrices and y stands for the Yukawa coupling. The 
scalar part of the action Sh in eqn. (|l|) is given by 

X,fl x 

+ £^ to-ll 2 , a = 1,2 . (4) 



The scalar fields <fi are complex, A is the scalar self coupling and k is the hopping parameter 
which is related to the mass parameter m c in the continuum formulation, m 2 c = (1 — 2A — 
8k) / k. The sums in eqs.@ and @ run over the lattice size A = L^L^.. 

The model defined in eqns.(|I]-[|) can be solved in the large fermion number (Nf) limit. 
It is convenient [0] to consider a modified form of the scalar action 



Sh = -K N J2 ( Pa(x)[(Pa(x + /j)+ip a (x-IJ,)}+J2 i Pa(x) 2 

x,n x 

+ X N (if 2 a (x) ~ N f ) 2 . (5) 
The usual lattice action eqns. (|l]-[|D is obtained by identifying 

k n = C 2 k , \ N = C 4 A , y N = Cy , (6) 
where the factor C satisfies the equation 

C 4 - (1 - 2X N N f )C 2 -2X N = (7) 
The scalar field (p a (x) is related to the original field (j) a {x) by 

<Pa{x) = <t>a{x)/C . (8) 
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In the 1/Nf expansion the couplings = yNfi/N, Xn = NfXjy are kept fixed to be 0(1). 
As Nf — ► oo the relations in eqn. (|6]) simplify, giving 



K N 



l-2\ N 



X 



N 



(i-2\ N y 



Vn 



1-2A 



X N < 



(9) 



N 



In the large Nf limit the constant mode of the scalar field dominates the path integral which 
suggests the Ansatz 



= JN f (a ,cp 2 (x)=0 



(10) 



where yNfa and yNfb correspond to the magnetization and staggered magnetization, re- 
spectively. The effective potential at leading order is 



— V eff (a, b) = -8K N (a 2 - b 2 ) + (a 2 + b 2 ) + X N {(a 2 + b 2 - l) 2 + Aa 2 b 2 ) - 2/(a, 6, 0) (11) 



where the zero temperature lattice integral is 

d 4 k 



/(a, 6,0) 



(27T)' 



log 



J^sin 2 kf, + y%(a 2 - b 2 



(12) 



2.2 Finite temperature 



A finite temperature on the lattice can be realized by a finite extension Lt in the temporal 
direction. The other directions on the lattice are kept large enough so that finite size effects 
are negligible. This leads to the usual finite temperature geometry of the lattice L\Lt with 
L s ^> Lt- The physical temperature is given in terms of the lattice spacing a 



T = l/L T a . 



(13) 



Varying the physical temperature can be achieved by changing Lt- Symmetry restoration on 
the lattice can then be detected in the following way: Assume that there is a phase transition 
from a symmetric to a broken phase at zero temperature. If now by heating the system, 
i.e. decreasing Lt, this phase transition shifts into the broken phase, we have found the 
symmetry restoration. Fixing the parameters of the theory, characterized by, for example, 
the renormalized scalar self coupling A r and the renormalized Yukawa coupling y r , to be 
in the broken phase at zero temperature, the system passes the finite temperature phase 
transition at some Lt corresponding to a critical physical temperature T c and for values of 
T > T c the system is in the symmetric phase. 

The lattice action of the finite temperature U(l) Higgs- Yukawa model is the same as in 
eqs.(|TJ-|p. The only change is that now the sums run over the finite temperature lattice with 
L s Lt- The steps of the large- Nf approximation are therefore very similar to the zero 
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temperature case. The only change is that at finite temperature the expression of I (a, b, 0), 
eq. (|T2|) , has to be replaced by 



/( fl) 6) L T ) = -L J2 [ 



d 3 k 



log 



/c 2 + 5>in 2 ki + f N (a 2 - b 2 



1,2,3 



(14) 



where &4 = 2-7r(n + 1)/Lt, n = 0, — 1. Note that in contrast to the continuum the sum 
over n is finite and can not be done explicitely. 



Replacing /(a, 6,0) by 7(a,6, T) in ([11]) gives the finite temperature effective potential. 
One may obtain different phases correponding to the locations of the minima of V e ff. 

(1) Symmetric (SYM) solution: There is a single minmum at a = b = 0. 

(2) Ferromagnetic (FM) solution: There is a minimum at a ^ 0, b = 

(3) Antiferromagnetic (AFM) solution: The minimum is at a = 0, b ^ 

(4) Ferrimagnetic (FI) solution: The minimum is at a ^ 0, b ^ 0. It can be shown that this 
solution does not exist for small values of the scalar self coupling A ~ 1. 

To determine the phase structure I have calculated the integral ( II ) which is needed for 
the effctive potential (0) numerically. I used mainly finite lattice sums as approximations 
to save computertime. However, I checked for various points that the lattice integrals give 
compatible results. 



As described in 11 



• I r mi j 



for large y values another type of large Nf expansion is possible. 
Here one keeps the Yukawa coupling ~ O(^JWf) and k n ~ 0(1/ J~N~f). As a result of this 
large- Nf expansion one obtains an effective action which is the AF-model in four dimensions 



S ef f = -K e// ^cr a (x) [a a (x + /i) + a a (x - /j,)) (15) 

1. The effective hopping parameter 



X,/l 



where the fields a a (x) have unit length, a a (x)a a (x 
K eff is given by 

Keff = K N Lp 2 + 



Nt 



(16) 



At finite temperature the steps of the calculation are the same with n e ff replaced by the 
one of the XY-model at finite temperature, K e ff(T). I determined k c (T) for the XY-model 
on lattices with Nt = 6 and Nt = 8 in time direction by a numerical simulation using the 
cluster algorithm. I find the phase transition still to be of second order and k c (T) shifted 
only slighly to larger k- values. This shift is less than 1% similar to case of the 0(4) model 
[H . Therefore the phase transition lines from the large Nf expansion in the strong Yukawa 
coupling region show only a tiny shift. 
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2.3 Numerical simulation 



I have performed numerical simulations at A = 0.0156, Nf = 2. Note that these parameters 
are the same as the ones used in Fig. 2a of ref. H] for zero temperature. The Hybrid Monte 



Carlo method |L2] was used for the dynamical fermion simulations. Each molecular dynamics 
trajectory consists of 10 steps with step size chosen such that the acceptance rate is around 
80%. As a check I did simulations at A = 0.0156 for Nf = 10 and found agreement with 
the large- iV/ expansion. To decide the order of the phase transition, I looked for hysteresis 
effects in the thermocycles. For each data point in the thermocycle about 50 trajectories are 
used as warmup and 100-200 trajectories are used in the measurement. 

The order parameters to detect the phasetransitions have been the magnetization v de- 
fined as 

V=<^ a >, ^a = J^Y,M^) (17) 

and the staggered magnetization v st given by 

Vst=<^ka>, ^=^£(-1)^^) , (18) 

where L is the linear size of the lattice. The measurements are done on 8 3 2 and a few on 
10 3 2 lattices. Although these lattices are certainly not sufficient to distinguish second order 
and weakly first order transitions, the combination and the agreement of the numerical and 
analytical results give a reliable determinination of the order of the phase transitions. 



3 Discussion and conclusion 



The results of the large Nf and Monte Carlo calculations are shown in fig.l for the Yukawa 
coupling y < 1.5. I have left out the results for the strong Yukawa region as they are 
indistinguishable from the zero temperature case. I checked explicitely also in this region 
that the Monte Carlo data agree with the large Nf predictions. 

In the figure the large-iVy results are shown as lines, where full lines indicate second 
order and dashed lines first order phase transitions. The Monte Carlo results for the phase 
transitions are exhibited as circles where open symbols mean first order and full symbols 
second order phase transitions. I find agreement between the Monte Carlo results and the 
theoretical prediction from the large Nf calculation not only for the position but also for the 
order of the phase transitions. There occur three different phases, a ferromagnetic (FM), a 
symmetric (SYM) and an antiferromagnetic (AFM) phase. As mentioned earlier, at large 
values of the Yukawa coupling a second symmetric phase appears, which is not shown in 
fig.l. The dotted lines in fig.l indicate the phase diagram at zero temperature as found in 
Although the structure and the order of the phase transitions remain the same at finite 
temperature, a clear shift to larger values of the Yukawa coupling is seen. 
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Let me concentrate on the phase transition between the ferromagnetic (FM) and the 
symmetric (SYM) phases which is relevant for the standard model. As discussed above, 
the observed shift to larger values of the Yukawa coupling indicates the expected symmetry 
restoration: A point (denoted as the star in fig. 1) chosen to be in the broken (FM) phase at 
zero temperature will be found in the symmetry restored phase at finite temperature. It is 
noteworthy that for the moderate temperatures that could be tested here the phasetransition 
is still second order and in the domain of attraction of the Gaussian fixed point. I conclude 
that symmetry restoration is not only an effect of a large N or pertubative apparoximation 
but survives also when non-perturbative methods are applied. 
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Fig.l The finite temperature phase diagram at A = 0.0156 and Nj = 2. The MC data 
are indicated by circles where the solid symbols denote second and the open symbols first 
order phase transitions. The solid and dashed lines are the results from the 1/Nf expansions, 
where the solid lines represent second order and the dashed line first order phase transitions. 
The phases are: (FM) ferromagnetic, (SYM) symmetric and (AFM) antiferromagnetic. The 
dotted lines indicate the zero temperature phase diagram || . Note that the phase transitions 
for finite temperature are shifted to larger values of the Yukawa coupling. This shows the 
expected finite temperature symmetry restoration as indicated by the starred point which is 
in the broken phase at zero and in the symmetric phase at finite temperature. 
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